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ABSTRACT 

Aims. Gamma ray bursts (GRBs) have great advantages for their huge burst energies, 
luminosities and high redshifts in probing the Universe. A few interesting luminosity 
correlations of GRBs have been used to test cosmology models. Especially, for a sub- 
sample of long GRBs with known redshifts and a plateau phase in the afterglow, a 
correlation between the end time of the plateau phase (in the GRB rest frame) and 
the corresponding X-ray luminosity has been found. 

Methods. In this paper, we re-analyze the subsample and found that a significantly 
tighter correlation exists when we add a third parameter, i.e. the isotropic 7-ray en- 
ergy release, into the consideration. We use the Markov chain Monte Carlo techniques 
to get the best-fit coefficients. 

Results. A new three-parameter correlation is found for the GRBs with an obvious 
plateau phase in the afterglow. The best fit correlation is found to be Lx oc T~'''^'' E?^'f^^. 
Additionally, both long and intermediate duration GRBs are consistent with the same 
three-parameter correlation equation. 

Conclusions. It is argued that the new three-parameter correlation is consistent with 
the hypothesis that the subsample of GRBs with a plateau phase in the afterglow be 
associated with the birth of rapidly rotating magnetars, and that the plateau be due to 
the continuous energy-injection from the magnetar. It is suggested that the newly born 
millisecond magnetars associated with GRBs might provide a good standard candle in 
the Universe. 
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1. Introduction 



Gamma-ray busts (GRBs) are one of the most powerful and energetic explosive events in 
the Universe. The observations of GRBs up to redshifts higher than 8 (Salvaterra et al. 2009; 
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Cucchiara et al. 2011) make GRBs to be among the farthest known astrophysical sources. 
Taking their considerable event rate into consideration, GRBs may be good candidates that 
can be used to probe our Universe. Several interesting correlations have been suggested for 
GRBs (Amati et al. 2002; Norris et al. 2000; Ghirlanda et al. 2004a; Liang & Zhang 2005; 
Dainotti et al. 2010; Qi & Lu 2010). Based on them, the cosmology parameters have been 
tentatively constrained (e.g., Fcnimore & Ramirez-Ruiz 2000; Schaefer 2003, 2007; Dai et 
al. 2004; Ghirlanda et al. 2004b, 2006; Amati et al. 2008; Wang & Dai 2006; Dainotti et 
al. 2008; Wang et al. 2009, 2011). 

To derive a meritorious constraint on the cosmology parameters, the most important 
thing is to find a credible standard candle relation for GRBs. Currently, no such a relation 
can be established when all GRBs are involved (Butler et al. 2009; Yu et al. 2009). The rea- 
son may be that different GRBs should be produced via various mechanisms. Interestingly, 
for a subsample of long GRBs with known redshifts and with a plateau phase in the after- 
glow, an anti-correlation has been reported to exist between the end time of the plateau 
phase (Ta, measured in the GRB rest frame) and the corresponding X-ray luminosity (Lx) 
at that moment (Dainotti et al. 2010, hereafter D2010). In this paper, we call Dainotti et 
al.'s two parameter correlation as the L-T correlation. The intrinsic scatter of this corre- 
lation is still too large to be directly applied as a redshift estimator (Dainotti et al. 2011). 
Additionally, normal long duration GRBs and the intermediate duration GRBs do not 
obey the same correlation equation (D2010), and the intermediate class seem to be more 
scattered in the plot. 

In this study, we have tried to add a third parameter, i.e. the isotropic 7-ray energy 
release (£J^,iso), into the correlation. We find that the new three-parameter correlation 
(designated as the L-T-E correlation) is much tighter than the previous L-T correlation. It 
is also obeyed by both the long GRBs and the intermediate calss. The L-T-E correlation 
may hopefully give a better measure for our Universe. In Section 2, we describe our GRB 
sample and the method of data analysis. Our results are presented in Section 3. Section 4 
is our discussion and conclusions. 

2. Sample & Data analysis 

According to Swift observations, many GRBs show a plateau phase in the early afterglow, 
prior to the normal power-law decay phase (Zhang et al. 2006; Nousek et al. 2006). In 
this study, we will mainly concentrate on the GRBs with such a characteristics. All our 
GRBs are taken from Dainotti et al.'s sample (D2010). In D2010's data table, totally 77 
GRBs are initially included, with known redshift and with a plateau phase in the afterglow 
light curve. After removing the intermediate class GRBs and some GRBs with relatively 
large errors, they finally limited their major statics to only 62 long GRBs. Here, we have 
re-selected the events by taking into account the following three criterions in our studies: 
(1) the plateau should be obvious (GRBs 050318, 050603, 060124, 060418, 061007, 070518 
and 071031 are removed by us, since their phatcau phase is not clear enough.); (2) the data 
in the plateau phase should be rich enough to show the profile of the plateau and the end 
time of the plateau as weh (GRBs 050820A, 060512, 060904 and 060124 are removed by us 
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due to this constraint.); and (3) there should be no flares during the plateau phase, since 
flares may affect the shape of the plateau light curve and lead to errors in the quantities 
that we are interested in (GRBs 050904, 050908, 060223A and 060526 are removed by us 

according to this condition.). As a result, our "golden sample" is consisted of 55 events 
in total, i.e., 47 long GRBs and 8 intermediate class GRBs (Intermediate class GRB is 
characterized by a short initial burst followed by an extended low intensity emission phase; 
Norris et al. 2006). The redshifts of our sample range from 0.08 to 8.26. 

For the end times of the plateau phase (Ta, in the GRB rest frame) and the X-ray 
afterglow luminosities at that moment (Lx = Lx(T'a)), we use the values of D2010. In 
D2010, Ta is derived through a phenomenological fltting model (Willingale et al. 2007), 
and Lx is derived from the following equation, 

47rDl{z)Fx 

= + 

where z is the redshift, Di,{z) is the luminosity distance. Fx is the observed flux by Swift— 
XRT at the end time of the plateau phase, and ^a is the spectral index of the X-ray 
afterglow (Evans et al. 2009). 

The isotropic 7-ray energy release in the prompt emission phase is 

£;^,iso = 47rD2(z)5boio/(l + z), (2) 

where 5boio is the bolometric flucncc, and can be taken from Wang et al. (2011). In the 
study of Wang ct al. (2011), ^boio is calculated from the observed energy spectrum $(£■) 
as (Schaefer 2007): 

where S is the observed fluence in units of erg ■ cm~^ for each GRB, and (i^min, i^max) 
are the detector threshold. The energy spectrum <^{E) is assumed to be the Band function 
(Band et al. 1993), 

\ BE^ otherwise 

where i?pcak is the peak energy of the spectrum, and a, [3 are the power-law indices for 
photon energies below or above the break energy respectively. At last, the complete data 
set of all our 55 GRBs are shown in Table 1, where the error bars are la range. 

We investigate if an intrinsic correlation exists between the three parameters of Lx, Tg, 
and E-y^iso as following, 

"'«<Ti?^)-« + '"og(^) + «Ml^). W 

where a, 6, and c are constants to be determined from the flt to the observational data. In 
this equation, a is the constant of the intercept, h and c are actually the power-law indices 
of time and energy when we approximate Lx as power-law functions of Ta and -E-^.iso- Due 
to the complexity of GRB sampling, an intrinsic scattering parameter, (Tint, is introduced 
in our analysis, as is usually done by other researchers (Reichart 2001; Guidorzi et al. 2006; 
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Amati et al 2008). This extra variable that follows a normal distribution of A''(0, af^^^) is 
engaged to represent all the contribution to Lx from other unknown hidden variables. 

To derive the best fit to the observational data with the above three-parameter correla- 
tion, we use the method presented in D'Agostini (2005). Here, for simplify, we first define 
xi = log(^), X2 = log( j|gj^), and y = log( io«CTg/s )- The joint likelihood function for 
the coeSicients of a, b, c and (Tint is (D'Agostini 2005) 

£(a,6,c,<Tint) oc JJ 



V r^r-nl {yi-a-bxi,i-cx2,if i 

where i is the corresponding serial number of GRBs in our sample. 

In order to get the best-fit coefficients, the so called Markov chain Monte Carlo tech- 
niques are used in our calculations. For each Markov chain, we generate 10^ samples ac- 
cording to the likelihood function. Then wc derive the the coefficients of a, b, c and CTint 
according to the statistical results of the samples. 

Our likelihood function can also be conveniently applied to the two-parameter L-T 
correlation case studied by D2010, by simply taking c = 0. We have checked our method by 
comparing our result for the L-T correlation with that of D2010. The results are generally 
consistent, which proves the reliability of our codes. 



3. Results 

In our study, we assume a fiat ACDM cosmology with H = 69.7 km • ■ Mpc~^ and 
Om = 0.291 (the same values as D2010). By using the method described in Section 2, we 
find that the best-fit correlation between Lx, Tg, and E^^iso is 

log(To^) = 1-17 - 0.87 log( j§5) + 0.88 log(^). (7) 

Figure 1 shows the above correlation. In this figure, the solid line is plotted from Eq. (7), 
and the points represent the 55 GRBs of our sample (the filled points correspond to the 47 
long GRBs and the hollow square points correspond to the 8 intermediate class GRBs). It 
is clearly shown that this three-parameter correlation is tight for all the 55 GRBs. 

Comparing Eqs. (5) and (7), wc find that the best values for the constants of a, b, 
and c in Eq. (5) are a — 1.17, b ~ —0.87, and c = 0.88 respectively. Figure 1 also clearly 
shows that there is still obvious scatter in the L-T-E correlation. To give a quantitative 
description of the scatter, we need to derive the la errors of these constants. 

The probability distributions of these constants as well as the intrinsic scattering pa- 
rameter ((Tint) are displayed in Figure 2. Prom this figure, we find that the probability 
distributions of these coefficients can be well fitted by Gauss functions. So we can easily 
get the la error bars for these parameters. Actually, the best values and the Icr errors for 
the coeflacients are a = 1.17±0.09, b = -0.87±0.09, c = 0.88±0.08, and a-.nt = 0.43±0.05, 
respectively. 

We have also explored the three-parameter correlation for all the 77 GRB events listed 
in D2010, using the same analytical method as for our "golden sample" of 55 GRBs. 
The best fit result is shown in Figure 3. The best parameter values and the Icr errors 
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Fig. 1. The best-fit correlation between ix, Tg, and -B-y^iso for our "golden sample". Y-axis 

is the X-ray luminosity at the end time of the plateau phase, i.e. Lx, in imits of 10^^ erg/s. 
Note that the X-axis is a combined quantity of Ta (in units of 10'^ s) and -B^.iso (in units of 
Xo53 erg), i.e. l.lT-O.STlogTa + O.SSlogi^^^^so. The filled points correspond to the observed 
data of 47 long GRBs and the hollow square points correspond to the 8 intermediate class 
GRBs. The solid line is plotted from Eq. (7), which is the best fit of the 55 observational 
data points. 

for the coefficients are a = 0.81 ± 0.07, b = -0.91 ± 0.09, c = 0.59 ± 0.05, and (Tint = 
1.15 ± 0.12. Comparing with the result of the "golden sample", although there is still an 
obvious correlation among Lx, and £'^,,130 for all the 77 GRBs, the intrinsic scatter of 
the L-T-E correlation is much larger now. However, it is very important to note that we 
exclude the 22 samples because they are most likely not physically belonging to the same 
group as the "golden sample" (for example, many of them do not have an obvious plateau 
phase), as judged from the three criterions in Section 2. 

In order to directly compare with the L-T correlation suggested by D2010, we have also 
fit the two-parameter correlation for our sample. The best-fit equation is 

log(To^) = (0-78 ± 0.14) - (1.16 ± 0.16) log( ^^j-). (8) 

This equation is consistent with the L-T correlation derived in D2010. Comparing Eq. (8) 
with Eq. (7) and from Figure 2, we find that the error bars of the constants in Eq. (8) (i.e. 
the L-T correlation) are generally significantly larger than those of Eq. (7) (i.e. the L-T-E 
correlation). Additionally, in the two-parameter fitting of Eq. (8), the intrinsic scatter is 
0.85 ±0.10, which is also markedly larger than that in the three-parameter correlation case 
(0.43 ±0.05). From the comparison, we see that the L-T-E correlation is really significantly 
tighter than the L-T correlation. 

For our GRB sample, we additionally find that the correlation coefficient of our L-T-E 
statistics is r = 0.92 and the chance probability is P = 1.05 x 10^^''. On the contrary, 
the correlation coefficient of the L-T statistics of the same sample is r = —0.73 and the 



6 



M. Xu & Y. F. Huang: L-T-E Correlation for GRBs with Plateau Phase 




-1.4 -1.2 -1.0 -0.8 -0.6 -0.4 0.2 0.3 0.4 G.b 0.6 0.7 0.8 

b 0:, 



Fig. 2. The probability distributions of the constants of a (upper left panel), b (upper 
right panel), c (lower left panel) in Eq. (5), and the probability distribution of the intrinsic 
scattering parameter Uint (lower right panel). According to these panels, the best values 
and the la errors for the coefficients are a = 1.17 ±0.09, b = -0.87 ±0.09, c = 0.88 ±0.08, 
o'int = 0.43 ± 0.05, respectively. 

corresponding chance probability is P = 5.55 x lO^'^. This also shows that the L-T-E 
correlation is much tighter than the L-T correlation. 

4. Discussion and Conclusions 

In this paper, a new three-parameter correlation is found for the GRBs with an obvious 
plateau phase in the afterglow. This L-T-E correlation is tighter than the L-T correlation 
reported in D2010. It has been shown that the intrinsic scattering of our L-T-E correlation 
is significantly smaller than that of the L-T correlation, and the correlation coefficient is 
correspondingly larger. However, we note that the intrinsic scatter of the L-T-E correlation 
is still larger than that of some correlations derived from prompt GRB emission (Guidorzi 
et al. 2006; Amati et al. 2008). In the future, more samples and more delicate selections 
might help to improve the result. 

The plateau phase (or the shallow decay segment) is an interesting characteristics of 
many GRB afterglows (Zhang et al. 2006; Nousek et al. 2006). This phenomenon can be 
explained as continuous energy injection from the central engine after the prompt burst 
(Rees & Meszaros 1998; Dai & Lu 1998; Zhang & Meszaros 2001; Dai 2004; Kobayashi & 
Zhang 2007; Yu & Dai 2007; Xu et al. 2009; Yu et al. 2010; Dah'Osso et al. 2011), or by the 
two component models (Corsi & Meszaros 2009), or by structured jets (Eichlcr & Granot 
2006; Granot et al. 2006; Panaitcscu 2007; Yamazaki 2009; Xu & Huang 2010), or even 
as due to dust scattering (Shao & Dai 2007; Shao et al. 2008). According to our L-T-E 
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Fig. 3. The best-fit correlation between Lx, and E^^i^o for all the 77 GRBs of D2010. 

The units of all physical quantity are the same as Figure 2. The X-axis is a combined 
quantity of 0.81 — 0.91 logTa + 0.59 log i?-y^iso. The filled points correspond to the observed 
data of 55 "golden" GRBs with error bars. The hollow diamonds correspond to 7 GRBs 
with too large error bars to be plotted in the figure, and the hollow circles correspond to 
other 15 discarded events. The solid line is the best fit for all the 77 data points. 

correlation (Eq. (7)), the X-ray luminosity at the end time of the plateau can be expressed 
as a function of the end time and the isotropic 7-ray energy release as, 

r ^ rp-0.87±0.09 77.0.88±0.08 /r,\ 
LxOCJa E^-^^ . (9j 

We believe that this relation can give useful constraint on the underlying physics. 

For the energy injection model, a natural mechanism is the dipole radiation from the 
spinning down of a magnctar at the center of the fireball. Note that the injected energy 
may not be Poynting flux, but can be electron-positron pairs (Dai 2004). These pairs 
interact with the fireball material, leading to the formation of a relativistic wind bubble. 
When the energy injection dominates the dynamical evolution of the external shock, the 
afterglow intensity should naturally be proportional to the energy injection power. So, Lx 
is actually a measure of the energy injection rate. According to Eq. (9), Lx is roughly 
inversely proportional to the timescale of the energy injection, Tg,. It hints that the energy 
reservoir should be roughly a constant. This is consistent with the energy injection model, 
which usually assumes that the central engine is a rapidly rotating millisecond magnetar. In 
different GRBs, the surface magnetic field intensities of the central magnetars may be quite 
different, leading to various energy injection luminosities and energy injection timescales. 
But the total energy available for energy injection is relatively constant (about rotational 
energy of the magnetar). It is mainly constrained by the limiting angular velocity of the 
magnctar, which again is determined by the equation of state of neutron stars. Additionally, 
according to Dai (2004), in order to produce an obvious plateau in the afterglow lightcurve, 
the total injected energy must be comparable to the original fireball energy (which may 
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be comparable to -E-y.iso)- This requirement is again roughly consistent with the item of 
-^°;fso'^°'°^ in Eq. (9). Based on the above analysises, we argued that the L-T-E correlation 
strongly supports the energy injection model of magnetars. It also indicates that the newly 
born millisecond magnetars associated with GRBs provide a good standard candle in our 
Universe. Thus the L-T-E correlation may potentially be used to test the cosmological 
models. 

Our sample contains 47 long GRBs and 8 intermediate class GRBs. From Figure 1, we 
see that both of these two classes are consistent with the same L-T-E correlation. Howerer, 
note that they behave very differently in frame work of the two-parameter L-T correlation. 
This is another important advantage of our three-parameter correlation. It indicates that 
magnetars may also form in intermediate class GRBs, and their limiting spinning is just 
similar to those magnetars born in long GRBs. A natural problem will be raised as to 
whether short GRBs with plateau phase in the afterglow also obey the same correlation. 
Unfortunately, the number of short GRBs meeting the requirement is currently too few. 

It is worth noting that many interesting physics could be involved in newly born mag- 
netars (Dall'Osso et al. 2009). The tops include the emission of gravitational waves, the 
cooling process, the evolution of the magnetic axis, etc. Some of the physics may affect the 
the energy injection process of the newly born magnetar delicately. We believe that further 
studies on the new three-parameter correlation may give useful constraints on the physics 
of newly born magnetars. 
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± 





10 


3.15 ± 0.1 


52 


79 


± 


0.013 


Long 


060502A 


1.51 


47 


27 


± 





19 


3.85 


± 


0.21 


52 


59 


± 


0.012 


IC 


060510B 


4.9 


47 


39 


± 





49 


3.78 


± 


0.48 


53 


64 


± 


0.011 


Long 


060522 


5.11 


48 


51 


± 





33 


2.07 


± 


0.31 


53 


05 


± 


0.026 


Long 


060604 


2.68 


47 


24 


± 





19 


3.98 


± 


0.18 


52 


21 


± 


0.069 


Long 


060605 


3.8 


47 


76 


± 





09 


3.48 


± 


0.08 


52 


66 


± 


0.034 


Long 


060607A 


3.08 


45 


68 


± 


2 


75 


4.14 


± 


0.02 


53 


12 


± 


0.012 


Long 


060614 


0.13 


43 


93 


± 





05 


5.01 


± 


0.05 


51 


32 


± 


0.006 


IC 


060707 


3.43 


48 


01 


± 





40 


2.94 


± 


0.36 


52 


93 


± 


0.025 


Long 


060714 


2.71 


48 


22 


± 





08 


3.11 


± 


0.07 


53 


06 


± 


0.016 


Long 


060729 


0.54 


46 


17 


± 





04 


4.73 


± 


0.04 


51 


69 


± 


0.021 


Long 


060814 


0.84 


46 


69 


± 





06 


4.01 


± 


0.06 


52 


97 


± 


0.004 


Long 


060906 


3.69 


47 


73 


± 





13 


3.62 


± 


0.12 


53 


26 


± 


0.042 


Long 


060908 


2.43 


48 


24 


± 





11 


2.46 


± 


0.09 


53 


03 


± 


0.010 


Long 


060912A 


0.94 


46 


37 


± 





23 


2.97 


± 


0.18 


51 


91 


± 


0.020 


IC 


061121 


1.31 


48 


35 


± 





10 


3 ± 0.09 


53 


47 


± 


0.004 


Long 


070110 


2.35 


48 


25 


± 





72 


1.89 ± 0.37 


52 


90 


± 


0.033 


Long 


070208 


1.17 


46 


88 


± 





15 


3.63 ± 0.14 


51 


58 


± 


0.060 


Long 


070306 


1.49 


47 


07 


± 





05 


4.42 ± 0.04 


53 


18 


± 


0.008 


Long 



Table 1. 55 GRBs of our sample. Data of the second, third, and forth columns are taken 
from D2010, the fifth column, i^-^jso, are calculated from Equation (2), where S'boio are 
taken from Wang et al. (2011). The last column is the type of the GRB, where Long means 
long GRB and IC is intermediate class GRB. All the error bars are la range. 
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(tIIB 






/( 










/('.-rg)] 


Type 


070506 


2.31 


47.63 


± 


1 


42 


2.87 ± 


1.42 


51.82 


± 


0.029 


Long 


070508 


0.82 


48.20 


± 





02 


2.75 ± 


0.02 


53.11 


± 


0.004 


Long 


070529 


2.5 


48.40 


± 





15 


2.34 ± 


0.15 


53.04 


± 


0.025 


Long 


070714B 


0.92 


46.85 


± 





20 


3.03 ± 


0.19 


52.30 


± 


0.033 


IC 


070721B 


3.63 


47.08 


± 





51 


3.58 ± 


0.51 


53.34 


± 


0.035 


Long 


070802 


2.45 


46.84 


± 


2 


72 


3.68 ± 


0.62 


51.96 


± 


0.047 


Long 


070809 


0.22 


44.15 


± 





76 


4.09 ± 


0.75 


49.43 


± 


0.062 


IC 


070810A 


2.17 


47.97 


± 





13 


2.83 ± 


0.12 


52.26 


± 


0.023 


IC 


071020 


2.15 


49.22 


± 





05 


1.84 ± 


0.05 


52.87 


± 


0.016 


Long 


080310 


2.42 


46.72 


± 





11 


4.08 ± 


0.11 


52.88 


± 


0.023 


Long 


080430 


0.77 


46.03 


± 





08 


4.29 ± 


0.08 


51.68 


± 


0.022 


Long 


080603B 


2.69 


48.88 


± 





26 


2.92 ± 


0.24 


53.07 


± 


0.011 


Long 


080810 


3.35 


48.24 


± 





08 


3.28 ± 


0.07 


53.42 


± 


0.031 


Long 


081008 


1.97 


47.79 


± 





24 


2.95 ± 


0.22 


52.85 


± 


0.047 


Long 


090423 


8.26 


48.48 


± 





11 


2.95 ± 0.1 


53.03 


± 


0.018 


Long 



